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Abstract. A linear Boltzmann equation is interpreted as the forward equation 
for the probabihty density of a Markov process {K{t),Y{t)) on (T x M), where 
T is the one-dimensional torus. K{t) is a autonomous reversible jump process, 
with waiting times between two jumps with finite expectation value but infinite 
variance. Y{t) is an additive functional of K, defined as v{K{s))ds, where 
\v\ ^ 1 for small k . We prove that the rescaled process N^^/'^Y{Nt) converge in 
distribution to a symmetric Levy process, stable with index a = 3/2. 



1. Introduction. 

The understanding of thermal conductance in both classical and quantum me- 
chanical systems is one of the fundamental problems of non-equilibrium statistical 
mechanics. A particular aspect that has attracted much interest is the observation 
that autonomous translation invariant systems in dimensions one and two exhibit 
anomalously large conductivity. The canonical example here is a chain of anhar- 
monic oscillators introduced by Fermi-Pasta-Ulam (FPU)[13J, for which numerical 
evidence shows a super-diffusive spreading of energy (see [IT9l1 for a general re- 
view). However, the rigorous analysis of energy transport mechanism presents 
serious mathematical difficulties and few results are obtained starting from micro- 
scopic dynamics. 

The canonical approach to this problem, starting with the pioneering work of 
Peierls [l23l] for the case of weak non-linearity, is to derive a Boltzmann-type equa- 
tion that will describe the energy transport in a kinetic limit. Recently, this ap- 
proach was carried out rigorously for weakly anharmonic FPU chains [i24i [H [2Tll . 
A linear Boltzmann equation was derived in [20] for the harmonic chain with 
random masses. The same linear Boltzmann equation appears also as limit of a 
random Schrodinger equation (see for example IfTTTl . IfTOH . Il25ll . Il2||). 

In [3] a kinetic limit was performed for a system of harmonic oscillators per- 
turbed by a conservative stochastic noise and the following linear Boltzmann equa- 
tion is deduced for the the energy density distribution of the normal modes, or 
phonons, characterized by a wave-number A; g[— 1/2, 1/2]: 

dtW{t,u,k) +v{k)duW{t,u,k) = [ dk' R{k,k')[W{t,u,k') - W{t,u,k)]. (1.1) 

The exact form of the scattering kernel R and the velocity v will be given below. 
The crucial features are, however, that the kernel R behaves like k'^ for small k, 
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'ik', and like k''^ for small k' , \/k, while \v{k)\ — > 1, as \k\ | 0. This conforms 
to the intuitive picture that phonons with wave number k travel with a velocity 
v{k) and are scattered with a rate R{k, k'). It is well known that super-diffusive 
spreading of energy is connected with the fact that the mean free path of phonons 
with small wave number k has a macroscopic length (ballistic transport), which 
follows essentially from the smallness of the rate with which these phonons are 
scattered, together with the fact that they travel with finite velocity. In fact, in [l3l] 
it is proved that this system exhibits anomalous conductance. 

To analyse the Boltzmann equation (II. ID is to exploit the fact that it can be 
reinterpreted as the forward equation for the probability density of a Markov pro- 
cess {K{t),Y{t)) on (T X M). Here K{t) is a reversible jump process with rate 
R and Y{t) is an additive functional of K, given as Y{t) = J^ds v{K{s)). In a 
phononic picture, the process Y{t) describes the trajectory of a phonon. For system 
with diffusive energy spreading, one expects that the law of the rescaled process 
Y {Nt) / \fN converges to the solution of heat equation. On the other hand, we 
expect that this is not true for systems exhibiting ballistic transport. 

To understand heuristically what is to be expected, it is convenient to introduce 
a discrete time Markov chain, Xi, that records the sequence of values assumed by 
the continuous time chain K{t) and the holding times r(Xi), i.e. the time the chain 
K{t) remains in the state visited in step number i. Then the process Y{t), at the 
time of the n-th jump of K{t), can be written as Sn = Ym=i 'T{Xi)v{Xi). We will 
see later that in our case, due to the fact that the transition kernel behaves as k'"^ 
for small k' with respect to the second argument, the stationary distribution of the 
chain Xi is of the form Tx{dk) = (j){k)dk, where (pik) ~ A;^ for A; << 1. On the other 
hand, the distribution of the holding time r(A;) is of the form P(r(A;) > s) ~ e~'^^^ . 
Hence, in the stationary distribution, we have that 

P(r(Xi) > s) ~ y 'K{dk)e-'^' ~ s'^/^ q_2) 

Thus, since v is antisymmetric and \v{Q)\ = 1, we expect r(Xj)f(Xj) to be in the 
domain of attraction of a stable law of index a = 3/2. Then, if the T{Xi)v{Xi) 
were independent random variables, n~^/'^Sn would converge to the corresponding 
stable law and n^^/°'S\.ns] to a stable Levy process. 

The corresponding problems of the convergence of dependent random variables 
with heavy tail distributions to a stable process has been studied extensively in the 
literature. A general, and very efficient, approach is explained e.g. in Durrett and 
Resnick [9]: first one uses methods from the theory of extremes of dependent ran- 
dom variables to study the convergence of the point process of scaled summands 
to a Poisson point process. Then one writes the sum as an integral with respect to 
the point process and uses a moment method to show that this integral converges 
to the corresponding integral with respect to the Poisson point process, which is a 
Levy process. 

The conditions needed to establish such a result are mainly required to assure the 
convergence to the Poisson process. Durrett and Resnick [9] express these in terms 
of the asymptotic behaviour of the conditional distribution. Davis [l6l] considers 
stationary processes and the corresponding well known mixing conditions (see 
[18]), while in [16] and [8] a non stationary generalisation is done, requiring 
some mixing conditions. Further results are to be found in e.g. [,7l [I4l [TSlI . 
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In this paper, we will follow the general strategy outlined above and use a cri- 
terion for Poisson convergence that can be found in in [5] and which has the 
advantage of being rather easily verified for an ergodic Markov chain. It requires 
asymptotic factorisation conditions for probabilities which hold on average. This 
will allow us to prove convergence first of the sum n^^/"S'[„t] to an 3/2-stable Levy 
process in the Ji-Skorokhod topology. Since Sn = Y{Tn), where T„ = XlILi '^(^j)' 
and since Y{t) is the piecewise linear interpolation of this function at the random 
sequence to times T„, once we will have shown that T[nt]/n converges to t, we will 
obtain that the rescaled process n^^/'^Y{ni) converges to the same limit process in 
the Skorokhod Mi -topology, which is the appropriate topology for the convergence 
of continuous process to a process with cadlag paths. 

Let us note that Jara et al [17] prove similar convergence results for additive 
functionals of Markov processes using Martingale and coupling methods and apply 
their results to the same model we consider here. Their methods do, however, only 
yield convergence of finite dimensional marginals. The methods we use here seem 
more straightforward and direct, and give stronger results. 

Finally, using convergence of n^^/°'Y{nt) to an a-stable Levy process, a G (1, 2), 
one can prove that the rescaled solution of the linear Boltzmann equation fll.lD 
converges to the solution of the following fractional diffusion equation 

dtW{t,u) = -{-dlr'^w{t,u). 

We refer to [flTll for the proof. Convergence of the rescaled solution of a linear 
Boltzmann equation to the solution of a fractional diffusion equation was indepen- 
dently proved by Mellet et al [.22,1 . with purely analytical techniques. 

2. The model 

We consider the process {K{t),Y{t)) described by equation (II. ID . Denoting by 
T the one-dimensional torus, we choose : T — > R and -R : T x T ^ as in [jSH , 
namely 

= (2.1) 

I sm(7rA;)| 

R[k, k') = ^ [2 sin^(27rA;) sin\nk') + 2 sin\2'Kk') sin'^ink) - sin'^{2'Kk) sin2(27rA;')] • 

(2.2) 

Observe that the rate kernel R is symmetric, not negative and it is equal to zero 
only if A; = or A;' = 0. We remark that despite the special case we consider, results 
depend essentially on the behaviour of v and M for small k, i.e. for A; -C 1 \v{k) \ ~ 1 
and R{k, k') ~ A;-^ VA;' G T/{0}, and R{k, k') ~ k'-"^ for A; < 1, VA; G T/{0}. The 
jump process K{t) is determined by the generator 

Lf{k) = [ dk' R{k, k')[f{k') - fik)], V/ G B{T). (2.3) 

The process Y{t) is an additive functional of K{t), defined as Y{t) = ds v{K{s)). 

Disregarding the time, the stochastic sequence {Xj}j>o of states visited by K{t) 
is a Markov chain with value in T, with a probability kernel P concentrated on T 
given by 

P{k, dk') = (j){k)-^R{k, k')dk' (2.4) 
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where </> : T ^ M+ is given by 

(Pik) = [ dk' R{k, k') = - sin2(7rA;)(l + 2 cos'^ink)). (2.5) 

We denote with P"^, m > 2, the m-th convolution integral of P. Since the prob- 
ability kernel P is regular and strictly positive and defined on a compact set, it 
is ergodic, i.e. there exists a strictly positive probability distribution tt such that 
V/c G T, P"(fc, ■) —>■ 7r(-), weakly, as n t oo. The stationary measure tt is given by 

7T{dk) = (f){k)dk. 

We define two functions of the Markov chain {Xj}j>o: the clock process, G 
M+j and the position process, Sn G M, by 

n-1 

T, = J2e^[HX,)]-\ (2.6) 

i=0 

and 

n-1 

Here {ej}j>o are i.i.d. exponential random variables with parameter 1. The clock 
process, T„, is the time of the n-the jump of the process K{t). It is a sum of positive 
random variables with finite expectation, as one can easily check using the explicit 
form of the probability density (see Eq. (15. 2D below) . The position process, Sn, is 
the value of the position of Y{t) at time T„, i.e. Sn = Y{Tn). It is a sum of real 
variables with zero mean and infinite variance. More precisely, for any i G N, for 
large A 

F\\e,<f){X,yMX^)\ > a] ~ X~'/^, (2.8) 
Let denote the right-continuous inverse of T„, i.e. let 

T-\t) = mi{n:Tn>t}. (2.9) 
We can represent the original processes, {K{t), Y{t)), as follows: 

Y{t) =S[T-^t)-l] + v{XiT-^t)-l]){t - TiT-i(t)-l])- 

In particular, Y{t) is the function defined by linear interpolation between its values 
Sn at the random points T„ (we take So = 0) . 

3. Main results. 

We assume that the initial distribution, jj,, of the process X satisfies the condition 

dfi{k)k~^ < oo, (3.1) 



which guarantees in particular that [eo0(Xo) ^] < oo. 
We define the rescaled processes 

T^(0) = ^V,j, 5a.(0) = ^V^j, T^\e) = ^T-\Ne), (3.2) 

where [-J denotes the lower integer part of -. Since r„ is a sum of positive vari- 
ables with finite expectation, we expect that both Tn{6) and T^^{6) converge in 
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probability (and thus in distribution) to 9, in the topology of uniform convergence 
on compact intervals. This will be proved in Proposition [S^ 

On the other hand, Sn is a sum of centred random variables whose tail behaviour 
is given in (I2.8D . Thus we expect that the rescaled process Sn converges to a stable 
process with index 3/2. This is the content of the following theorem. 

Theorem 3.1. Let Sn he the process defined in fl3.2l) . Then for any < T < oo, the 
process {SN{6)}o<e<T converges to {V{6)}Q<0<r, where V is a symmetric Levy process 
stable with index 3/2. Convergence is in distribution on the Skorokhod space ofcddlag 
functions equipped with the Ji — topology. 

Combining this theorem with Proposition 18.21 we will prove that Sn ° T^^ con- 
verges in distribution to V . This will imply our main theorem. 

Theorem 3.2. Let Sn, T^^ be the processes defined in (\3.2^ . For every < T < oo, 
the process {SN{T^^{6))}Q<e<T, converges to {V{6)}Q<Q<r, where V is a symmetric 
Levy process stable with index 3/2. Convergence is in distribution on the Skorokhod 
space of cadldg functions equipped with the Ji — topology. 

Moreover, for every N > 1 the process YN{t) = ds v{Ks) is the function 

defined by linear interpolation between its values Sn{0) at points Tn{0), with 6 G 
[0, T^^(T)]. In particular, Yn converges to V in distribution in the Mi-Skorokhod 
topology. 

4. Sketch of the proof of Theorem 13.11 

In this section we present the key steps of the proof of Theorem 1. The technical 
details will be given in Sections 5, 6, and 7. As we mentioned in the introduction, 
we follow the strategy of considering the sequence of the point processes associated 
to Sn- At the first step we define this sequence of point processes and we show 
that it converge to a Poisson point process. Then we prove that the limit process 
for Sn exists and it is a Levy process stable with index 3/2. Finally, we prove the 
tightness for the sequence Sn- 

4.1. Point processes. Define the real valued random variables 

= (piXnrMXn), U G Nq. (4.1) 

For some fixed c > 0, we decompose SNifi) into two parts, SNid) = S^{6) + 
5<(^), where 

[Ne]-i 

"=° (4.2) 



|<cAr2/3}. 



n=0 



We will see later that S§{6) vanishes as oo and then c — > 0. More precisely: 
Lemma 4.1. Let S^{9), 9 G [0, T] be the process defined in (\4.2\ l. Then for every 

E\\S§{e)\^]<Coe^/d + Ci9N-^/^, (4.3) 
where Co, Ci < oo are positive constant. 
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We will prove this lemma in the next section. 

On the other hand, will be connected to two Poisson processes. 

We split the sum S^{6) into two parts: 

n=0 

+ 7^ l{e„^„<-ciV2/3}. (4.4) 



n=0 



Defining the random variables 



^i,N = j;pj^ l{v.>o}, Xi,N = j;pj^ l{^,<o}, (4.5) 
with values in M+, M_, respectively, and the associated point processes 7^^ and 

we can rewrite S^{6) as 

POO rd p—c 

S^{9) = / / xnjj{ds,dx) + / / xn^{ds,dx) (4.7) 

Jo Jo J-oo 

The following Proposition states that the two point processes TZ^ and TZ]^ converge 
to Poisson point processes. 

Proposition 4.2. Let TZ^, TZ~^ he the point processes defined in (f4.6P . Then the point 
process TZ% converges in distribution to Poisson point processes TZ^ on M+ x W^, resp. 
M+ X M_, with intensity measures dt x di'^{x) = ^a\x\~^^'^dx, with a > some explicit 
constant. 

We will prove this proposition in Section 6. 
4.2. Limit process for Sn- We define a process V on IR+ with values in M by 

pB pod p9 pQ 

V{e)= / / TZ+{ds,dx)x+ / / n-{ds,dx)x, (4.8) 



-'0 -'0 ^-oo 



and 

f9 /-oo 



/•y POO py r~'^ 

V>{e)= j j n^{ds,dx)x+ j j n-{ds,dx)x, (4.9) 

Jo Jc Jo J -oo 

where 7Z^, 7^ are the Poisson point processes defined in Proposition 14. 2i 

Proposition 4.3. Let he V, the process defined in (\4.8h V is well defined and it is 
symmetric a Levy process stahle of index 3/ 2, with Laplace functional \E'(A) = Co|Ap/^, 
where Cq = a J dy\y\~^/'^{l — cosy). 

Moreover, Vc ^ V as c I 0, in the topology of uniform convergence on compact 
intervals. 
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Proof. We split V{e) into two parts V{0) = V>{0) + V<{6). Clearly, the two pro- 
cesses, if they exist, are independent. Since, for any compact interval, /, the total 
intensity of / x (c, oo) equals 

/peso 
ds J du^{x) = \I\ac-^/^ < oo (4.10) 

the process (6) is a finite sum, almost surely. 

Moreover, by direct computation and using the fact that E[Vc'^(6')] = 0, we get 

p8 PC I'd pO 

E[\V<{e)\^]= ds du+{x)x^+ ds du~{x)x^ = 69ac^/\ (4.11) 

Jo Jo Jo J-c 

hence also V'^{6) is almost surely finite, and tends to zero in probability, as c | 0. 

Then V{6) is almost surely finite, and since it is the sum of two independent pro- 
cesses which are right-continuous and have independents increments, it satisfies 
the hypothesis of a Levy process, and it is full characterized by the one-dimensional 
distribution. This is uniquely determined by its characteristic exponent \I/ : M — > C 
defined as 

g-t*(A) _ ^ [exp{iXV{9)}] (4.12) 
which by direct computation is given by 

/oo 
diy{x) (1 - e^^^) , (4.13) 
-oo 

with the Levy measure z/ defined on M \ {0} given by du^x) = a\x\~^^^dx. For every 
XeR 

vP(A) = \Xf/'a [ dy \yr'''\l - cosy) = Co\X\'/', (4.14) 
Jr 

which is the characteristic exponent of a s)aTimetric Levy process stable with index 
3/2. □ 

Proposition 14.31 together with Lemma 14.11 and Proposition 14.21 implies conver- 
gence of the finite dimensional distributions of Sn to V. Since V has stationary 
and independent increments, it is in fact enough to prove the convergence of the 
one-dimensional distributions. 

Corollary 4.4. For any 9 e M+, 

SN{e) ^ V{9), (4.15) 

in distribution. 

Proof. From the representation fl4.6D and Proposition 14. 2[ and since the intensity 
of [0, 6] X (c, oo) is finite, it follows readily that 

SNid) ^ K(^), as iV T OO, (4.16) 

in distribution. Moreover, Vc{9) ^ V{9), as c I 0, and 

limlimsup|^^(^)| =0, (4.17) 

in probability, by Lemma \4A\ This implies the assertion of the corollary. □ 

To conclude the proof of Theorem 1, we need to complement this corollary with 
the proof of tightness of the sequence Sn- This will be postponed to Section 6. 
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5. Moment estimations. 

In this section we collect and prove some useful moment estimates. 
We start with some preliminary results on the transition probability density for 
the Markov chain {X„}. They are given by 

p{k,k') = (f){k)-^R{k,k') (5.1) 

= 8 \; sin^ ttA;' + 8 , , sin^ vrA;' cos^ vrA;' • 

1 + 2cos2(7rA;) ^ ^ 1 + 2 cos2(7rA;) ^ ^ ^ ^ 

We denote hy p"^{k, k') the m-step transition densities, i.e. 

p\k,k')=p{k,-)op{.,k') = [ dk,p{k,k,)p{h,k') 

Jt 

and in the same way, for every m> 1 

p^^\k, k') = p^{k, ■) o p{., k') = p{k, ■) o k'). 

Observe that P"^{k, dk') = p'^{k, k')dk' , where is the m-th convolution integral 
of the probability kernel P. 

In the next proposition we give an explicit formula for p™. 

Proposition 5.1. For every m > 1, p'^{k,k') has the following form 

pm^k, k') =8—^^p^^^^ [a^ sin^rck') + sin^rck') cos^nk' 
1 + 2 cos^(7rA;) 

Sill (ttJv^ r A. 1 /\ 1 2/ 7 /\ 2/ 7 / 

+8:r— t; Tr~r\ sm {nk ) + dm sm (vrfc ) cos [nk 

1 + 2 cos^{7Tkj '- 

where 

ai = di = 1, 6i = ci = 0, m = 1 

<^ 0^171 7 7 } dfn ^ 1 , Tfl ^ 2 

Proof By direct computation, 

p\k, k') =8^^^M.— U sm\7ik') + 62 sin2(7rfc') cos^T^k' 
1 + 2 cos^ [-Kk) '- 

sm r 7T r 4/7/\ 7 '2, / 1 f \ ( t I 

+8:; r,, , x \c2 sin (vrfc ) + ct2 sin (vrfc ) cos (vrfc 

1 + 2 cos^(7rA;) '- 



(5.2) 



(5.3) 



(5.4) 



m > 2 (5.5) 



with 02, &2, C2, d2 positive, a2 = d^. In the same way for m > 2 we find expression 
(15. 2D . where the coefficients a^, 6^, c^, are given by the following recursive 
formula 

Cm+1 <^m+l / \ C2 d2 J \ Cm 

In particular a^, &m, Cm, are positive for every m > 2 . Moreover, by condition 

J.J, d/c' p^ik, k') = 1, Vm > 1, we get the relations 

which says that the coefficients are uniformly bounded, i.e. a^, bm, Cm, dm < 1 for 
every m > 1 . □ 

We now give the proof of Lemma 14. 1[ 
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Proof. (Lemma [4. ID . Let us write 



E 



[Af6»J-l 



n=0 



[Ne]-i 



(5.6) 



{e„|V„|<cAr2/3}l{e„|^„|<cAr2/3}J 



n=l 0<m<n 

Let us focus on the first sum on the right hand side. For every n> Iwe have 



E 



I I ^ 

\en'^n\ l{e„|i/;„|<cAf2/3} 



where 



dze l_dfi{ko) j__dk p'\kQ,k)\'ilj{k)\ l{,\^(k)\<cNV^}^ 

3| cos(7rfc)| 



\^{k)\ = \<p{kr'vm 



2sin2(7rA;)(l + 2 cos2(7rA;)) ' 



(5.7) 



(5.8) 



Using the exphcit form of given in 05. 2D . an elementary computation reveals 
that 



E 



(5.9) 



with Co, Ci < oo are positive constants. Thus 



E |e„^„| l|,„|^„|<,^2/3} <Coe^c + Ci9N'^'\ (5.10) 



n=0 



The second sum on the r.h.s. of Eq. (15. 6D is in fact equal to zero. Namely, for all 
n > m> 1 we have 



E [en'^Pnem^ipm'^{e„\■4,„\<cN■2/3}'^{e,r^\1pm\<cN^^'■^}] 



dz ze 



duue-"" / dfi{ko) / dkp"'{ko,k)^{k)l 



(5.11) 

{\i>(k)\<cN^/3z-^} 



X / dkp^"" '")(A;,fc')^(^')l{|V'(fc')l<c^^''^«- 



□ 



since tp{k) = -ipi-k) and p\k, k') = p\k, -k') \/k, k' eT,\/i> 1. This concludes 
the proof of Lemma [4TT1 

The following related lemma will be needed in Section 8. 

Lemma 5.2. There exists C < oo, such that for n > 1 

Proof For every n > 1 we have 



(5.12) 



dz ze ^ dfiiko) / c^/c p"(A;o, ^) | l{2|v>(fc)|>Af2/3}, 

(5.13) 

and using the explicit formula (15. 2D for one easily finds that this expression is 
bounded by CiV-^/^, with C < oo. □ 
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Lemma 5.3. Let S^{9), 9 e [0,T], be the process defined in (\7.12\) . Then for every 

< r < s <t <T 



E [\S^{s) - S<{r)\'\S<{t) - S<{s)\'\ < Ao{t - r)^ 
for some constant Aq < oo. 
Proof For every 0<r<s<t<T we have: 

\SUs) - SUr)\'\SUt) - S<{s)\' 

^ lNs\~l lNs\~l 

i=[Nr\ j=|_AfrJ 
[AftJ-1 [AftJ-1 

^ ehV'hl{eh|^^|<Ar2/3}e;V'/l{e,|^,|<Af2/3}- 
h=lNs\ l=lNs\ 

By (15. IIP , the only terms with non-vanishing expectation are 

[ATsJ-l lNt\-l 



(5.14) 



(5.15) 



i=lNr\ 

For every h > i we have 



|<Ar2/3} ^ {eh'iphl l{ehl^/'hl<Ar2/3}- 
h= [Ns] 



E 



leil/ji] l{e.|^-|<Ar2/3}|e/iV^/i| l{eh|^/)h|<iV2/3} 



(5.16) 



(5.17) 



dz z e 



2„-z 



duu e 



2-u 



dfi{ko) 



X / dkp\ko,k)\^{k)\H{,\^^k)\<NV^} / dk' p^'' ^)(A;,/c')|V'(^')riH^(fe')l<^2/n- 
Using (Q, we find that for all A; G T, 

rfA;>('^-)(A;,A;')|V^(fc')plH^(,0l<iv2/3} (5.18) 



'u 



cos^i^nk) sin^(7rA;) 
1 + 2cos2(7rA;) ^ 1 + 2 cos2(7rA;) 

where B is a finite constant. Thus, Wko G T 

dk p\ko, A;)|V'(A;)pl{^|^(fc)|<;v2/3} / dk' p^^-'\k, A;')lV^(^')Pl{«|,^(fc')l<Jv2/3} 



< 



\ LatS/s^ / cos^(7rfco) sin^(7rfco) \ 

fz y/u \ 1 + 2 cos^(7rA;o) 1 + 2 cos^(7rA;o) / 



(5.19) 



with i?i < oo. Finally we get 



E 



\S<{s) - S<{r)\'\S<{t) - S<{s)[ 



< 



1 



< Ao{t-r) 



i=lNr\ h=lNs\ 
2 



which is the assertion of the lemma. 



(5.20) 

□ 
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6. Proof of Proposition 14.21 

This section is devoted to verify conditions of the Theorem 2.1 in QS]], which 
guarantees the convergence of the point processes defined in (14. 6D to Poisson point 
processes. We recall the statement of Theorem 2.1 of [5]. 

Theorem 6.1. [[5]] Denote by Xla^vW ordered sequences of 

different indices with values in {1,..,A^}. Let be {Z^^^} (resp. {Z^j^}), 

N > 1, 1 < i < N, an array of random variables with values in (resp. in 
Assume that for every £ > and all sets of positive constants (ri, .., r^) 



J2 P [\ztN\ > n, 



I7± I 



3/2 ■■■ 3/2 



as 



oo, 



(6.1) 



for some a > 0. Then the point process J2iLo^ ^i/N,z^j^ 



converges in distribution to a 
Poisson point process 7^^ on [0,T] x M± with intensity measure dt x dv^{x), where 



dv^ix) 



\a\x\ 



X. 



Our goal is to verify these conditions for the random variables {X^j^}, {X-j^} 
defined in 04. 5D . In order to do so, we need the following lemma. 

Lemma 6.2. For every i > 0, for every sequence of different indices, ii,..,ie G 
{1, .., A^}^ and all sets of positive constants, (ri, .., t^), the following statements hold: 

(i) if {ij+i — ij) > 2 for every j = !,..,£— 1, then 



(6.2) 



with Ci > 0. 
(ii) If (ij+i - ij) = I for some j = 1,.., 



1, then 



lim N'T^^\..Tf^^ F[eM> N^/'n,...,eMJ > iV^/V,] = 0. 

N^oo 



fe^V2 j,i2 
Proof Let us consider 



dz e"" / dh p'^ik,k')l 



{|V(fc')|>Af2/^r2-i}! 



(6.3) 



(6.4) 



with m > 1. Using (I5.2D and the explicit expression for ijj 05. 8D , one easily finds 
that for m > 2 this quantity is bounded by 



cos^(7rA;) 



+ 



l + 2cos2(7rA;) 1 + 2 cos2(7rA;) 
with Ci < cx). On the other hand, for m = 1, we get 



1 1 



(6.5) 



poo p 

j dz j dk' p{k,k')l^\^(^^.,^\^j^2/3, 
Jo Jt 



'0 

<Co 



cos (tt/c) 



+ 



sin^(7rA;) 



l + 2cos2(7rA;) 1 + 2 cos2(7r/c) 



(6.6) 



with Co < CO. 
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Next we consider 

due'"' j c^fci p''(A;,/i;i)l{i^(fc^)|>^2/3^^„-i} / dz e 

X / rfA;>'"(/ci,A;')l||^(fc/)|>^2/3^2^-i|. 
By (1631) . (I6^ we find that for n,m>2 



(6.7) 



due P"(fc,fci)l{|^(fci)|>w2/3^i«-i} / 



(6.8) 



cos^(7rA;) ^ sin^(7rA;) 



1 + 2cos2(7rA;) 1 + 2 cos2(7rA;) 



1 1 1 

jV2 3/2 3/2 ^ 



with C2 > 0, while if m = 1 and/or n = 1, then (I6.7D is of order By repeating 

this procedure, one finds that for every ni, .., > 2 the following inequality holds: 



Jj Jo 

p poo 

X / c^A;2P"'(A;i,A;2)l{|^(fc2)|>iv2/3r2.2-i} / c^^rge" 



^2 



T Jo 



dz£e~ 



(6.9) 



X / dhp'''{h.i,ke)l^^^^f^^^^^^2/Zr,z-^} 



- iV^""^ 3/2 ••• 3/2 



1 + 2cos2(7rA;) 1 + 2 cos2(7rfc)_ 



where Ci are finite constants for ni, .., > 2. If = 1 for some i, then the 1. h. s. 
of (Ia9l) is of order o(iV-^). 
Assume without loss of generality, ii < 12 < .. < ie.h Then 

P POO p 

l^dfi{ko) I dzie-'' j^dh p''{ko, h)!^^^^^^^-^^^ ^2/ 3^^,-1 y 

(6.10) 



T 
00 



X / dzo e ''■'^ 



T Jo 



dk^v"^ *'(^i,^2)l{i^(fc2)i>^2/3,^,-i| / dz^e 

J 



□ 



The proof of the lemma is now just an application of (16. 9D . 

Recalling the definition of Xf^, X^^ given in G4.5D , by the symmetry of the 
probability density we have 



(6.11) 



Let us denote by J2f3pf{e) ^^"^ ordered sequences of different indices 

[h, ..,ie] with values in {1, ..,N}. Then Eaj,{e) = ^' Ep^iey 
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We choose ii < 12 < ... < ie and by denoting with mi = ii, rrij = ij — ij^i, 
Wj =2,..,i, we get 



N (6.12) 

mi , . . — 1 
7n-i--{-..~\'m£<N 



Fixed M <^ N, we spHt the sum on mi into two parts: 



N 

, . . ,m£ — 1 

M N 

l+...+mi)\ > 

(6.13) 

mi=l m2,..,m^^l 

N-e+i N 

m-|^ + ..+m£<iV 



Using Lemma [621 we find that the first sum on the r. h. s. of (I6.13D is bounded by 
CeT^'^^"^ ...T^^^'^M/N, with Ce > 0. Let us consider the second sum on the r. h. s. of 
(I6.13D . We spUt the sum on m2 into two parts and we get 



Af-^+l N 



^ ^ P [eil^^J > N^/^Tu...,ei\^^m,+...-,m,)\ > iV'/V,] (6.14) 



mi=M + l m2...,m^-l 

N-e+l N-i+l N 



Y Yl Y ^ [eM^m,\ > iV2/Vi,...,Q|^(^,+...+„,o| > ^'^'^^] 

=A/+1 m2=l 

+0{M/N). 



mi=M+l m2=M-\-l m3,..,m^-l 
m,-\-\-..-\-mf<N 



By repeating this procedure for all the sums, finally we get 



N 

J2 P hk-il > N^^''n,.:,ei\4Jim,+...+m,)\ > iV^/V,] (6.15) 

m]^ + . .-\-m£<N 

N 

m ]^ , . . , m£ — A / + 1 
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Now we show that in the last expression we can replace the probability P with 
the invariant measure vr. We have 



N 



mi ,. . ,m£ — M-\- 1 



(6.16) 



^/t Jo Jt 



Let us consider the second sum on the rh.s. of (I6.16D . We have 

OO f 

dze / dke b"''(^^-b - l{|^(fc,)l>7v2/3r,.-i} 



(6.17) 



< sup {<Pikr'\p"''{k,k')-<p{k')\} / dze / c?fc,0(A;,)l{|^(,^)|>^./3_-i| 



k,k'€T 



which is finite since p'^{k, k')(j){k') ^ G C°°(T x T) for all m > 1. Moreover, since 
by direct computation 



OO 

dz 



dk 0(fc)l{|^(A;)|>Ar2/3^^-i} < Cq — ^^, 



with Co > 0, using Lemma (621 we find that the second sum in (I6.16D is bounded, 
for N large enough, by 

^^~m^^' c^4j^-4n ^^v {m-'\p"'ik,k')-<p{k')\} 

(6.18) 

<Ci^...^ sup snp {<P{kr'\p"'{k,k')-<j>{k')\}. 

T^' Tp m>M+l fc,fc'GT 

Now we consider the first sum of (16.16ft . By repeating this procedure £-times, we 
can replace the transition probability densities p™(A;, k'), i = 1, — 1, with the 
invariant density (l){k'). This gives an error S^ii, M,ti, ..,Ti), which satisfies the 
inequality 

\S^{i,M,T,,..,n)\<Ce4j^...4j^ sup sup {<P{k')-'\p'^{k,k')-<l){k')\}. 

T^' Tp m>M+l fc.fc'eT 

(6.19) 

By ergodicity, for every i>l and all sets of constants r,; > 0, i = 1, .., £ 

Mm \SN{i,M,Ti,..,Te)\=0, ViV G N. (6.20) 

M—*oo 
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Then, for large enough. 



J2 P hl^^il > iV'/Vi,...,Q|v^(„,+...+„,)| > N^/^n] (6.21) 



m-^-\-. .-\-m£ <N 

N 



Y] ^i=i( [ dz /c/fc0(A;)l{i^(fc)|>;v2/3r,^-i} ) 



^1 1 ■ 



+^,v(^,M,ri,..,r,) + 0(M/iV) 



— - L + 0(iV-5/3) 



AT^ 3/2 ••• 3/2 

^1 



+ ^^(£,M,ri,..,r,) + 0(M/iV). 



We take the Hmit A^, M oo with MjN 0. The proof of the proposition follows 
by (I6.20D and by relation 

N 

lim V N-^ = i/e\ 



m-^ , . . ,m£ — l 
m2^'{-..'{-m£<N 



7. Tightness. 

In this section we conclude the proof of Theorem 13.11 by proving tightness of 
the sequence S^. This is relatively easy due to the strong convergence properties 
stemming from the weak convergence of the point processes 71^. 

As criterion for tightness in the Ji topology we use slight variant of Theorem 
13.2 (with condition 13.5 replaced by 13.8) from [H]. 

We define the modulo of continuity on D 

Wf{6) = sup{(|/(t) - /(ti)| A \f{h) - fm ■.h<t<t2,h-h< S}. (7.1) 

The sequence {Pn} of probability measures on {D, V) is tight in the Ji -topology if 
and only if 

(i) For each positive e there exist r such that 

Pnlf ■■\\f\\oo>r]<e, n>l. (7.2) 

(ii) For each e > and r] > 0, there exist 6 > 0, and a integer no such that 

Vn > no 

Pn[f : Wf{6) >v]<e, n> no, (7.3) 

and 

Pn[f : Wf{6) >v]<e, (7.4) 

Pn[f:\fi6)-f{0)\>v]<e, (7.5) 

Pn[f:\f{T)-f{T-6)\>r]]<e. (7.6) 

In order to verify these conditions, we start with some definitions and prelimi- 
nary results. For every < r < s < t < T, let us denote by 

nirst = min (|5'Ar(s) - SN{r)\, |S'Ar(t) - 5'Ar(s)|) (7.7) 

and 

Lt = sup rrirst (7.8) 

0<r<s<t<T 
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The following inequalities holds (see H, (10.4), (10.6)): 



sup \Sn 

eg[o,T] 



sup \Sn 

ee[o,T] 



<Lr+\SN{T)\, 

<3Lr+ max N^'^^^eMi 

0<i<[NT]-l 



Fix > 1. We have 



P 

where 



sup Sn{0) > T 
ee[o,T] 



< P 



sup S^{9) >t/2 
ee[o,T] 



P 



sup S^{9) > r/2 

ee[o,T] 



i=0 
1 '^^'"^ 



i=0 



Let us consider inequality (I7.9D with Sn replaced by S^. We have 



P 



sup S^{e) > r/2 

9g[0,T] 



< P [L^ > r/4] + P \S^iT)\ > t/4: . 

By direct computation, for all < r < s < t < T we have by Lemma 15.31 

PK<, > r/4] < P \\S<{s) - S<{r)\\SUt) - ^^(^)| > r'/W 



A 



Then we can use Theorem 10.3 of Q to get 



B 



P [L^ > r/4] < -T 



(7.9) 
(7.10) 



(7.11) 



(7.12) 



(7.13) 



(7.14) 



(7.15) 



where i? is a constant. Moreover, by inequality (I4.3D with c = 1, 6 = T we get 



P 



\S<{T)\>t/4 



16^ 
< -E 



(7.16) 



with some constant D. Finally by fl7.13D , fl7.15D and (I7.16D , we get that for every 

> 1 



P 



In order to estimate P 



sup 

9e[o,r] 



sup S^{e) > r/2 
ee[o,r] 



B^2 



(7.17) 



sup6ig[o.T] 'S'^(^) > , we observe that 

1 t^^'"^ 



i=0 
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Then by a straightforward computation (see Lemma [5^2]) 

INT\-1 



P 



sup S^{e) > r/2 
ee[o,T] 



2 1 
< r- 

- T iV2/3 



J2 E[e#i|l 



{e.|V>l>JV2/3}J 



i=0 



(7.18) 



- T- jV2/3 jVl/3 



2K.T- 



with constants Kq . This inequaHty with G7.11I) and (I7.17D prove (I7.2D . 
Now we prove (I7.4D . Again, for each > 



P [^5^(5) > r/] < P Ws<{S) >r,/2 + P ^^>(5) > r,/2 



(7.19) 



To estimate P 



Wg<{5) >r]/2 , we divide [0,T] in intervals of length 5 > 0. If 
u^5<(5) > ^/2, then exists k, 1 < k < T /6, such that 

sup jmin (|5<(s) - 5<(r)|, |^<(t) - 5<(s)|) , (A; - 1)5 < r < s < t < (A: + l)^} 
> V2. 



Observe that 



Thus, by (f7J4D . 



min ( [^^(s) - ^^(r)!, \S'^{t) - S'^{s)\) = m, 



rst' 



1 



P[m<,> V2] < A-(t-r) 



Hence 



(A;-l)5 <r <t < (A; + 1)5. 



Co 



(7.20) 



(7.21) 



P [sup {m< t : (A; - 1)5 < r < s < t < (A; + 1)5} > r]/2] < -^6% 
with some constant Co (see Theorem 10.3 in H). Using this result we get 

¥[ws<{5)>v/2 

[r/s] 

< ^ P [sup {m<j : (A; - 1)5 < r < s < t < (A; + 1)5} > t]/2] (7.22) 



fc=i 



<^r5 



with some constant Ci. 

Now let us consider the quantity 



rrirst = min \S^{s) - 5jv(r) , \S^{t) - S^{s) 
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for < r < s < t < T, with t - r < 5. If m>t > 0, then eM > A^^/^ ej\^pj\ > N'^/^ 
for some [rN] < i < [sN] - 1, [sN] <j< [tN] - 1. We get 



P 



Ws>i6)>r]/2 



<¥[3i,j < NT, \i-j\< 25N : al^^i] > N'^^^ ej\ipj\ > N"^^^] 

[Nr]~l i+[2NS] 
i=0 j=i+l 

< CT5 



(7.23) 



where the last inequahty follows by Lemma 16.21 Equation (I7.23D together with 
(17:221) proves (17^41) . 

The proofs of equations Q7.5D , (17.6D are easier. We give only the proof of (17. 6D , 
since the other is similar. We have 



F[\SN{T)-SN{T-6)\>r]] 



< P 



S<{T)-S<{T-S)\>v/2 +P \S>{r)-S>{T-6)\>v/2 



(7.24) 



where, using a second moment estimate (essentially the same proof of Lemma 
14TT]) , we get 



Moreover, 



P 



\s<{r)-s<{r-s)\>v/2 

INT\-1 



P 



1 



J2 e„V^nl{|e„^„ 



|<Af2/3} 



n=lN{r-5)\ 



> W2 



4 1 ^ 
< E 



[ATTJ-l 
n=[Af(r-5)J 



2n 



P 
< P 



< 



S>{T)-S>{r-S)\>v/2 

^ INT]-1 

r^TJ E l{len^„|>iV2/3} > ^/2 

n=[Af(r-5)J 

[7vrj-i 



Ar2/3 
2 1 



n=[Ar(r-<5)J 



(7.25) 



(7.26) 



and this together with (17.25!) proves (17. 5D . 
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8. Proof Theorem [3T2] 

Let us consider the two rescaled processes Tn, T^^ defined in (I3.2D . We want to 
prove the convergence in probability of both processes to the function 6. We start 
with the following Lemma. 



Lemma 8.1. Let be the process defined in l\3.2\) . Then Ve > 

lim P \\TN(e) -e\>e]=0, G [0, T] 



(8.1) 



Proof. Le us denote by the expectation value with respect to the invariant mea- 
sure TT. Since [ej</)(Xj)"^] = = 0, 1, then 



P [\TN{e)-e\ >e]<F 



[Af6»J-l 



^1 J2 (e.0M-i-E. [e,0(X,)-^]) 

i=0 



> e 



N 



We denote hy e' = e — N ^, which is positive for large enough. Let us introduce 
the following notations: 



= ei(j){Xi) ^l|g^^(jY^-)-i<^2/3}, rr^ = ei(f){Xi) ^l{g^^(jY^)-i>^2/3}, 



(8.2) 



Vi = 0, 1,.... Wehave 



P 



[Af6»J-l 



1| J2 {e^<P{X,)-'-E^[e,^{X,)-']) 



< P 



i=0 



> 6' 



N 



(r< - V 



i=0 



> 



+ P 



[7V6»J-1 



i=0 



> 



(8.3) 



By a first moment estimation (see Lemma [5721) . we get 

[JVSJ -1 



E 



TV- 



i=0 



1 



iVl/3' 



with < oo, thus we can neglect the second term on the rh.s. of (18. 3D . For the 
first term we have 



P 



< 



N 
4 1 



[Nei-i 

i=0 

lNe\-i 



- y 

Ml ^ 



< 11 



> 



^1/3n 



+ 



4 1 



\Ne\-\ 



E[r,<r<-E.[T,<]E.[r<]] 



(8.4) 



£2 iV2 

where, using a second moment estimate (see for example (I5.10D ). 



4 1^ 



i=0 



(8.5) 
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Let US consider the second sum in 08. 4D . For fixed 1 < M < A^, we split it into 
three parts: 

[Af6»J-l 

[Ndl^l M-1 [Ne\-1 

= E E[r<r<-E.[r<]E4r<]] + J] E [r<r< - [r<] E4r,<]] 



ij=0 
l<|j-i|<M 



1=0 j=o 

\j-i\>M 



[Ne]-1 [Af6»J-l 



i=M J=o 

b-i\>M 



(8.6) 



For every i,m> 1, we have 



E [^^^+m] = / dz ze ' duue dn{ko) / dk p\ko, k)(j){k) ^l{^<^(fc)-i<iv2/3 



X / dk' p"'{k,k')(j){k') ^l{,^(fcO-l<iV2/3}. 



Using (15. 2D one finds that E [r^r^^] ~ thus the first and the second sum in 
(|8^ are 0{M/N). Moreover since Tcidk) = (f){k)dk, one get 

= E^[t<t<^^] + / rf2;e"^z / d/i(A;o) / dk p'{ko, k)(p{ky^l{,^(^kyi^N2/3y 
Jo Jt Jt 

poo f 

X / du e-^u / dk' (p^ik, k') - 0(A;')) (l)ik')-\uHk')-^<Ny^ 
Jo Jt 

+ / dze-'z / c?/i(A;o) / dk {p\ko, k) - (j){k)) ^{ky^l^^^^^yi^^i/^yE^ [t<^^] , 
Jo Jt Jt 

and then Vi, m > 1 

E [r<T;< J - E. [r<] E^ [7;< J < C /" rffcmax sup {0(fc)->'(fco, fc) - • 
Thus the third sum in (I8.6D is bounded by 

^-2^2^ / dkmaxsnp {(l){k)-^p\ko,k) - ^{k)\] , 
Jt *>a^ fcoeT 

which, by ergodicity, goes to zero for M ^ oo. We prove the Lemma choosing 
N,M ^oo, with M/N -^0. □ 

Now we prove that the processes Tn{9), Tj^^{9) converge in probabihty (and thus 
in distribution) to the function 9. 

Proposition 8.2. Let T^, T^^ he the processes defined in dOD . Then Ve > 0, V5 > 
3iVo such that WN > Nq 



sup {\TN{e)-0\}>e 

6»6[0,r] 



<6, 



(8.7) 
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and 



sup {\T^\e) - e\} > e 



< s. 



(8.8) 



Proof. We give the proof only for (18. 8D . the proof of fl8.7D is similar. 

Fixed rj > 0, we divide [0, T] in [T/rj] +1 intervals of length less or equal than 
rj. We observe that 



P 



sup {\T^\e) - e\} > e 

0e[o,T] 



< P 



sup {T^\e)-e} >e 

9e[0,T] 



+ P 



sup {e-T^\e)}>e 

eelo,T] 



(8.9) 



We have 

P 



sup {T-\9)-9} > 
0e[o,T] 

ir/vi 



ir/v] 

k=0 



sup {T-\e)-e}>e 



(8.10) 



< 5^ P [{T^'lik + - (A; + l)r^} > e - v] 

k=0 

= ^ P [T^\kr]) -kr]> e-T]] 

k=l 

where we used the fact that T^^ is non-decreasing. For every e > 0, we choose rj 
such that e' = e — ?7 is positive, and we rewrite the relation T^^ikr]) — kt] > e' as 

T-\Nk7]) > Nkr] + Ne', (8.11) 

where T^^{t) is defined in (I2.9D . By definition the following relation holds: 

T[T~^{Nkr))~l] < Nkrj < T[T--L(Nkr))], (8.12) 

then inequahty (18. IIP imphes 

Nkri > TiNkri+Ne'-i\, (8.13) 

and thus 

J2 P [Tj,\kr^) - kr^ > e'] 

k=l 

< ^ P [{ki] + e' - N-^) - TN{kr] + e' - A^"^) > e' - A^"^] . 

k=l 



(8.14) 



In the same way one can easily prove that 



P 



sup {9-T^\9)}>e 

6»e[0,T] 



ir/r,\ 



< ^ P [TN{kr] - e') - {kr] - e') > e'] . (8.15) 



fc=0 



By Lemma [8II1 Ve' > 0, W > 0, 3No such that F \\Tn{0) - ^| > e'] < S', ViV > No, 
for all 9 e [0,T], and this conclude the proof of (fCT . □ 
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In order to prove that Sn ° T^^ converges in distribution to the Levy process 
V, we have just to show that the distance between the two process Sn, Sn ° T^^ 
goes to zero in probabiHty (see Theorem 3.1 in O). By denoting with p(-, ■) the 
distance in the Skorokhod Ji topology. 



lim F[p{SN,SNoT^')>e] <P 

TV — >oo 



sup |S'Ar(AAr(t)) 

te[o,r] 



SNoT^\t)\>e 



(8.16) 



where \n 



[0, T] [0, T] is a sequence of increasing homeomorphisms such that 
^ 0. We have 



P 



sup \SN{\N(t)) - Sn oT^'^(t)\ > e 



< P 



sup |S^(A^(t))-^^or/(t)| >eln<^ sup \T^\t)-t\<(j 

te[o,T] te[o,T] 



+ P 



sup |r^^ 



(t) -t\> a 



for every a > 0. The second term on the r.h.s. goes to zero for ^ oo (see 
Proposition [8^, while the first term is bounded by 

P[sup{|5jv(Ajv(t)) -5jv(s)|,|t-s| <a, t,se [0,T]} > e] , 

which goes to zero for the tightness (see G7.4D ) and the fact that HAa? — / ||oo ^0. 

It remains to prove the convergence of the process Yn to the stable process V. 
The basic idea is that the step-function sequence {SN{6)}o<e<r and the continuous- 
time sequence, given by the linear interpolation of {Sn{6)}, are asymptotically 
equivalent, i.e. if either converges in distribution as N — oo, then so should 
the other, and they should have the same limit. This is proved, for example, in 
Il26ll . Section 6.2. Then one can easily extend this result to {5'Af(T^^)}o<6)<r and 
{yN}o<e<T, since the waiting time between two jumps has finite average. 
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